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Abstract 

The capillary instability of a LC cylinder in magnetic field is considered using an energy approach. 
The boundary problem is solved in the linear approximation of the anisotropy Xa of the magnetic 
susceptibility x- The effect of anisotropy, in the region 1 ^> \x\ > \Xa\ ^S> x 2 , can be strong enough 
to counteract and even reverse the tendency of the field to enhance stabilization by enlarging the 
cut-off wave number k s beyond the conventional one set by Rayleigh. 
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1 Introduction 



Theoretical predictions of the continuum theory of the nematic liquid crystals (LC) were successfully 
confirmed in many experimental observations [1]. One of the most studied effect is the influence of an 
external field on the orientational distribution of the LC director field n(r). The physics involved in 
a competition between the mechanical and field forces that can be developed in LC media, calls for 
careful analysis. Well known example is the Freedericksz effect [1], observed when nematic LC cell with 
initial uniform distribution of n(r) is subjected to an external magnetic field. In many cases, the stabi- 
lizing elastic forces compete with the destabilizing magnetic field giving rise to a critical phenomenon. 
However, the critical phenomena in LC can be sustained even if both the elastic and magnetic fields are 
defined as stabilizing. Such critical phenomena exist due to the effect of surface tension which tends to 
minimize the surface area by distorting the initial shape of the system. 

The nematicity of LC's, being a source for elastic properties, results in enhancement of stability 
of LC jets [2], as compared to ordinary liquids. A similar stability enhancement appears in ordinary 
liquid jets with isotropic magnetic permeability (see [3] and [4]) when they are subjected to an external 
magnetic field. Unlike the elasticity, the external field has a critical value beyond which instability of the 
jet is completely suppressed for all disturbance wavelengths [4]. Nematic LC's are usually anisotropic 
diamagnetics with positive anisotropy Xa of the magnetic susceptibility [1]. This poses an additional 
challenge with respect to the above mentioned phenomena. It is reflected by the extra terms in the LC 
hydrodynamics due to orientational interaction between the magnetic field H and LC director n. 

The static version of capillary instability in liquid jets is known as the Plateau instability in the 
liquid cylinders [2], [4]. It dates back to the classical works of J. Plateau [5] who defined the problem of 
finding a surface of liquid with a minimal area S given its boundary d£l at fixed volume. The problem 
relates to the principle of minimum free energy at equilibrium. Further generalization is called for if the 
excess free energy W, of the cylinder, comprises different types of energy that reflect a more complex 
structure of the liquid (e.g. elasticity [2] etc), as well as its capacity to interact with external fields. 

The purpose of this work is to extend the theory of Plateau instability in LC cylinders [2] so as to 
include the effect of static magnetic fields. Here the motivation is both experimental and theoretical. 
Experimentally, the question is how to set the initial orientation of director n collinear with the LC 
cylinder axis. A weak magnetic field can serve to this end. Theoretically, the framework outlined in [2] 
can be extended so as to incorporate the influence of external fields on the evolution and stability of the 
LC cylinder. In particular, the Plateau instability is studied with respect to the effect of the magnetic 
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anisotropy of the LC cylinder. 



2 Free energy of LC cylinder in the presence of magnetic field 

Consider an isothermal incompressible LC cylinder in a uniform magnetic field Ho that is applied in 
free space along the cylinder axis. We assume a rigid boundary condition (BC) where the director is 
tangentially anchored at the free surface of the LC cylinder. The magnetic susceptibility tensor \ of the 
LC is assumed anisotropic, symmetric, and independent of the magnetic field. In the reference frame 
related to the cylinder axis its diagonal terms are X||>X_i_ while its off-diagonal term is Xan- Z n r , where 
ft 2, f fly* 0X6 the axial and radial components of director n, respectively. 



Figure 1: Undisturbed (left) and disturbed (right) homotropic LC cylinder subjected to an axial and 
uniform external magnetic field. 

When the LC cylinder (assumed long compared to its diameter) is undisturbed the total free energy F° 
of the system is defined by 

cyl 

where the integral represents the volume itR 2 L of the undisturbed cylinder which is enclosed by the 
surface dQ® yl . The term £® = o~ Jqqo ds = 2iraRL stands for the surface free energy of the undis- 
turbed cylinder, where R, L and a denote its radius, length and surface tension respectively, and no 
is permeability of a free space. We specify the commonly used harmonic deformation of the cylinder 
as r(z) = R + Co cos fez, where k = 2ir/A, A being the disturbance wavelength. Let the extent of 
deformation be characterized by a length (o, such that (o/R = e <^ 1. 

Deformation of the cylinder shape changes the magnetic field H(r) over all space M 3 , while the 
director field n(r) is changed only within the internal domain fl cy i- Following Plateau, we assume 
conservation of the cylinder volume 
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The total free energy F of the disturbed cylinder takes the following form, 

F = £ 8 + S n + S% + E% . (3) 

The first term in (3), which stands for the interfacial energy of the disturbed cylinder, is classically 
known due to Plateau [5] 

r nC 2 
£ s = a ds = 2iraRL + <r-%L (k 2 R 2 - l) . (4) 

Jdn cyl 2# 

The second term in (3) is due to the elastic deformation of the director field n(r), and in the single 
elastic approximation is given by 

K f 

£ n = — (div 2 n + rot 2 n) dv , (5) 

2 Jn cyl 

where K is the elastic modulus. The last two terms in (3) correspond to the effect of the magnetic fields 
in the internal Q cy i and external 1R 3 \ Q C yl domains 

= IT / H > ~ T f '"^ £ 5 = ! rf H > - T I 6X H 2 (r)d,, (6) 



L cyl ~ Jii cy l ~ JIK"\il cyl - JIK"\il cy l 

where m H(r) and ex H(r) are the internal and external magnetic fields, respectively. If the deviations of 
the director n = n° + n 1 from its initial orientation n° along the z direction are small, then 

n° = 0, n° = l, 1 > n\ > \n\\ ~ {nlf . (7) 

The magnetic energy density in the second term of £^ ( scaled by hq/2) reads 

fi jk ' m Hj' m H k = {l + X \\){' m H z ) 2 + (l + x±)C"H r ) 2 + 2 X an 1 r n° z '"H r i "H z , (8) 

where Hj k is the LC relative permeability tensor: [i zz = 1 + x\\, Hrr = 1 + X±, A*zr = Xa^r n z- The 
excess free energy W of the system is defined as, 

W = F — F° . (9) 

From the mathematical standpoint, the variational problem for minimization of W, supplemented with 
constraint (2) for all smooth surfaces dft cy i, is known as the isoperimetric problem. The cylinder 
instability can be studied assuming small perturbation in its shape. In this case the Plateau problem 
becomes solvable in closed form. The fields m H(r) and ex H(r), which must satisfy Maxwell equations, 
can be presented as small perturbations of Ho, 

in H(r) = H + in H 1 (r) = (H + m H l z , ' m H l r ) , ex H(r) = H + ex H 1 (r) = (H + ex #i, ex H}) , (10) 



where according to the assumption e <C 1 the following approximations apply 

{ m H}., "Hi, ]n Hl, "Hi} = {'X, e X, ''"hi, e *hl} x cH . (11) 

The dimensionless fields m,ex /i^ z (r, z) are dependent on the coordinates as indicated By virtue of trans- 
lational invariance of the problem 

in %(^ + A)= in ' ex /4>^) (12) 

we set L = A and evaluate the free energy per unit wave length. Substituting (1), (4)-(8) and (10) into 
(9) we obtain in the e 2 -approximation 

V = t^-oR (k 2 R 2 - 1) + £ f (div 2 n + rot 2 n) dv - ^ U , (13) 
L 2 2L J Qcy[ 2L 

where the magnetic part was calculated in Appendix A 

U = f Ul + Xll )(i»Hl) 2 + (l + x± )(>"Hl) 2 }dv+ [ {("Hl) 2 + ("Hl) 2 \dv + 
Jn cyl <- ' JR3\n cyl L J 

2H [ X a[ n r n z ' n Hldv + (l + X \\) [ in Hldv+[ ex tf>). (14) 



3 Boundary problem and its solution 

The magnetostatics of the disturbed LC cylinder is governed by Maxwell equations for the internal 
in H(r) and external ex H(r) magnetic fields and the Euler-Lagrange equation apply for the director field 
n(r), 

rot in H = rot ex H = , div in B = div ex B = , (15) 

{ dr d(d r n r ) dz d(d z n r ) dn r } ^ n ~^ ^ ' ^ x dx ^ ^ 

where m B and ex B denote internal and external magnetic inductions, respectively 

m Bj = fionjk m H k , fi jk = (1 + x±) $jk + Xarijn k , ex Bj = fi ex Hj , Xa = X\\ ~ X± • ( 17 ) 

Equations (15), (16) must be supplemented with boundary conditions (BC) at the interface r = R, 

( in H,t) = rH,t), ( in B,e) = ( ex B,e), (e, n) = , (18) 

where t and e stand for tangential and normal unit vectors to the surface, respectively. Since the surface 
deformation is small, linearization can be applied, 

t r = -e z = dC/dz , t z = e r = ^l-(d(/dz) 2 ~1. (19) 



A standard way to solve the problem is to introduce a director potential 0(r) and two magnetic poten- 
tials 3>in(r) and 3>ex(r) as follows 

nj = — , '"H 1 ^) = -V$ in , ex H 1 (r) = -V$ex, |V$ in |, |V$ex| < flo . (20) 

$in(r) and <£ e x(r) satisfy the first two equations in (15). The last two equations in (15) yield, 

d 2 <S>- Id / <9G\ <9 2 <3> 

(1 + X||) ^ + (1 + X.) A 2 <l> in = Xa H Q \l (r£) , ^ + A 2 * ex = , (21) 

where A 2 = d 2 /dr 2 + 1/rd/dr is the two-dimensional Laplacian. The variational equation (16) gives, 
Making use of the commutation rules 



we finally arrive at 



A - — — - — (a — ^ 9 

2 r 2 <9z 2 / <9r <9r y 2 <9z 2 / ' 



d 2 \^ ~ /• ,9 2 ^ _, , 



A 2 + — j <D ex = , (^A 2 + q 2 ^-2 1 $ in = (a 2 - l)H A 2 e , (23) 

9z 2 y k i + x ± 

BC (18) can be reformulated as follows 

dz dz ' dr dr dz ' <9r dz 

A weak decoupling of the equations (23) makes it possible to solve the boundary problem in closed 
form. Assuming 

{<S>- m (r,z),$ ex (r,z),e{r,z)} = {(j); n (r), 4> ex (r), 9(r)} xsinkz (25) 

we find 

(A 2 - k 2 ) 0ex = , (A 2 - a 2 k 2 ) </> in = Xa H A 2 e , (A 2 - k 2 ) 9 = 2^>|% in , (26) 
The following BC exist at r = R 

^ex = 0in, (1+X±)^-^=X^C0, f r =~kC0. (27) 

The two last equations in (26) can be represented through the determinant equation 

[A 2 - D 1 A 2 + D ] ^ ^ j = ^ j , where (28) 

D X = (2 + Xa )fe 2 + Xa5 2 . D = (1 + Xa )fc 4 , g 2 = . (29) 



Factorization of the differential operator in (28) gives 

A 2 - L>iA 2 + D = (A 2 - If) (A 2 - if) , where \\ 2 = -(p 1 ± y/vj , (30) 
V = Dj -4D = Xa[xa(k 4 + g 4 ) + 2(2 + X a)k 2 g 2 ] . (31) 

The fundamental solutions of (28) which are finite at r = are the following 

Mr) = cj^hr) + c 2 I (l 2 r), 9{r) = bj^hr) + 6 2 / (Z 2 r), ^^-bj = lj ~ " * Cj -, j = 1, 2, (32) 

where I m (x) is a modified Bessel function of the 1st kind and order m, and Cj,bj are indeterminates. 
The first equation in (26) produces a solution in the exterior domain which is finite at r = oo 

(f>ex(r) = c 3 K (kr) , (33) 

where K m (x) is a modified Bessel function of the 2nd kind and order m. All indeterminates ci,c 2 ,C3 
can be found from BC (27) by substitution therein the expressions (32), (33). Further simplification 
comes after substitution of C3 = c\Iq(1\R) / 'Ko(kR) + c 2 Iq(1 2 R) I 'Ko(kR) and using 1 + x± — 1 



ci 



hh{hR) + k%$$I (hRj\ + C2 [| 2 / 1 (z 2 i2) + fc^gg/ (z 2 12)l = x± Ko^o 



K (kR)- 



(34) 



where the identities I (x) = h(x) and K' (x) = —K\{x) for the derivatives were used. Straightforward 
calculations give 

r V ■ l 2 — a 2 k 2 

Cj = k( H -±, b j = k( ^ J yt- , J = 1,2, (35) 

1 1 XJj 

where Fj,j = 0, 1, 2 are determinants of (2 x 2) matrices 

r = a (I 2 - I 2 ) h{hR)h{l 2 R) + [(ah - l 2 )I (hR)h(l 2 R) - (al 2 - h)h(l 2 R)h(hR)} 

-rr = (x^-xWk^^ + x.k^^Io^R), 

r 2 = (xA-xW^'-^ + xak^^HhR). (36) 

3.1 Xa^expansion of the solutions 

The complexity of expressions (32) in conjunction with (35) and (36) makes further evaluation of the 
problem excessively difficult. Therefore, we develop in this Section another approach for solution of 
the amplitude equations (26) endowed with BC (27). Bearing in mind that for most nematic LCs the 
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anisotropy is small \x a \ < |x±|> |X||I 1 we seek the linear in Xa representation of the functions <f>\ n {r) and 
9{r) 

Mr) = ^i„(r) + XaMr) , 9{r) = 9(r) + Xa 9(r) . (37) 
The isotropic parts 4>- m (r), 9(r) together with the external potential </> ex , satisfy the following equations 

(A 2 - k 2 ) ex = (A 2 - k 2 ) ^ in = (A 2 - k 2 ) 9 = , (38) 

supllemented with the BC at r = R 

<Aex = in , (l + X±)^-^ = X±H k(o, ^ = -K0. (39) 

The solutions 4>\„(r), 9(r), <j>ex(r) of (38) were found in [2] and [4] 

^ in (r) = Aiio(fcr), 0(r) = A 2 I (kr) , c/> ex (r) = A 3 K (kr) , (40) 

where the coefficients C« are given by, 

A 1 = CoX±H kR ^ R) , A 2 = --^- A 3 = (ox±H kR *°% R) , (41) 
T(kR,x±) h{kR) T(kR,x±) 

and T(x, a) = l+axIi(x)Ko(x). Henceforth, T(kR, X ±) — 1 in accordance with xj_ <C 1. The amplitude 

equations for the remaining functions 4>\r,(r) and 9{r) can be found by inserting (37) into (26) and making 

use of (38) 

(A 2 - k 2 ) in = Vo(fcr) . ^ = {M + A 2 ^ )fc 2 = "|0y . (42) 
(A 2 - k 2 ) 9 = A e h{kr) , A e = 2A 1 ^ = 2 x± C kRK (kR)f^- . 
The BC for in (r) and 6{r) is, 

MR) = 6(R)=0. (43) 
After simple calculations (see Appendix B) we obtain 

Mkr) = ^±G(kr)I (kr) , 6{kr) = ^G{kr)I Q {kr) , (44) 

where 



x The characteristic magnitudes of the magnetic susceptibility % an d its anisotropy x« f° r the classical nematic LC's 
4-metoxybenziliden-4-butilanilin (MBBA) and para-azoxyanisole (PAA) can be found in [6] : x± ~ X\\ ~ — 10 _s and 



Xa « icr 6 . 
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Recalling the definition (20) of potentials <I>j n (r), ( I ) ex( r ) an d O(r) we get 



nl(r,z) = kA 2 h(kr) ll + Xa 



A e 
A 2 k 2 
At 
A x k 2 
A, 



G(kr) + 



G(kr) + 



G' r (kr) I (kr) 
k h{kr) 

G'(kr)I (kr)~ 
k Ii(kr) 



sin kz 
sinkz . 



- m H l z (r,z) = kA 1 I (kr) jl + Xa -^G(kr) j cos kz , 
ex H^(r, z) = kA 3 K 1 (kr) sin kz , ex ^(r, z) = -kA 3 K (kr) cos kz . 



where 



A e „ /i H 2 R 2 h(kR)K Q (kR) A a 

Xa-TTo - ~ 2 XaX± T7 7~n > Xa~ 



1 



(46) 

(47) 

(48) 
(49) 

(50) 



A 2 k 2 ' ' WUL K kR ' A\k 2 x±kRh(kR)K Q (kR) ' 

Before going to further calculation of the excess free energy W by (13), let us estimate the main 
contributions of the anisotropy Xa to the distribution of the fields n(r, z) and m H 1 (r, z) in accordance 
with (49) and (50). 



First, as follows from Figure 2, G{kr) and G(kr) + 



G'(kr) I (kr) 
k h(kr) 



which are both continuous monotone 




Figure 2: Plots of the functions G(kr) (solid curve) and G(kr) + G ^ j"^) (dashed curve) 



growing funstions, are bounded as follows 



-0.21 < G(kr) < , and 0.75 < G(kr) + < OS) . 



(51) 



k I\(kr 

In order to simplify further calculation we consider, henceforth, these functions as constant N\ and N 2 
respectively 



r(h s AT rfh s , G'(kr) J (fer) 
G(fer) = iVi , G(fcr) H r — - - — = N 2 , 



(52) 



k I\(kr) 

The next simplification comes for nl(r,z). Indeed, bearing in mind (50) we conclude that for the 
influence of Xa on n\.(r, z) to be significant a huge magnetic field Hq is required 



H > H. 



1 K 



V\XaX±\ R V W 



(53) 



The magnitude of this field can be as high as 10 8 A/m for classical LC's with radius R ps 10/xm. In 
fields which are significantly lower than H,, the behaviour of the director n(r) is dictated primarily by 
competition between bulk elasticity and surface tension of the LC's, and goverened by dimensionless 
parameter x= K/aR [2]. Recasting (46), (47) and (48) gives 

nl(r,z) = kA 2 h(kr) sin kz , ex H^(r, z) = kA 3 Ki(kr) sin kz , ex #i(r, z) = -kA 3 K (kr) cos kz , 
- m H}(r,z) = kA 1 I 1 (kr) (l + Xa -^N^J smkz , - in #i(r, z) = kA 1 I (kr) (l + Xa -^N^J cos kz. 

Now we are in position to calculate the magnetic part of the excess free energy W according to (14), in 
the limit Xa < X± ^ 1 ( see Appendix A). 

\2_. u^oikR) ( _ 1T)J , un ^ T s f1 .n\ , oAr 



-U = KH(oH yx±kRj^ \x±kRh(kR)K (kR) + 2iVi^J • (54) 
The elastic part of W was found in [2] 

£ n = itLKk 2 d ■ ( 55 ) 

Inserting (54) and (55) into (13) we get 

W = ^(^) 2 -f(kR, X ,H ), (56) 

where 

f(kR, x,H ) = (kRf(l + 2x) - 1 + ^ g °™ fc ( ^ ) {xlkRh(kR)K (kR) + 2AT lXa ) . (57) 

All the terms in (57), except the last one, describe the stabilization of a LC cylinder due to the existence 
of isotropic susceptibility irrespective of its sign [4] and due to the elasticity of the LC phase [2]. The 
influence of the last term in (57) which accounts for the anisotropy of x can De significant and even 
dominating. The latter occurs if 

\Xa\ » xl , (58) 

as indeed is the case in classical LC materials (MBBA, PA A). Here the physical situation changes 
completely. The cylinder is destabilized with the corresponding cut-off k s R 

k s R^\ + 2 X a\N 1 \^^. (59) 
a 

The interesting property of (59) is the fact that the cut-off k s R extends beyond the range (0 < k s R < 1) 
of the classical Rayleigh instability. This kind of extension cannot be obtained as a field or elastic effects 
in the absence of anisotropy Xa of the magnetic susceptibility. 
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4 Conclusion 



• The capillary instability of a LC cylinder in magnetic field is considered using an energy approach. 
The excess free energy, which includes terms due to surface, LC's elasticity, and magnetic field 
is used to find extremum conditions associated with instability. The boundary problem is solved 
and then expanded in terms of the anisotropy Xa of the magnetic susceptibility. 

• The excess magnetic free energy, which was founded to be a function of the isotropic susceptibility 
(%) squared proved to have an anisotropic part linear in Xa- This indicates that the effect of 
anisotropy can turn dominant provided that 1 S> > \Xa\ ^ X±- This means that the effect 
of anisotropy can be strong enough to counteract and even reverse the tendency of the field to 
enhance stabilization by extending the cut-off k s wave number beyond the conventional range set 
by Rayleigh. 

• As the existence of magnetic anistropy is not limited to complex fluids such as LC's, the result of 
this work can be considered of a more general nature. 
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A Contribution of the magnetic field to free energy 

Evaluate the contribution hqU /2 of the magnetic field inside Q cy i and outside M 3 \ £l cy i of the disturbed 
liquid cylinder to the excess free energy W 



U = 



dv + 



(Al) 



2 X 



[ ( ] "H z ) 2 dv- [ H 2 dv + j U>"H z ) 2 + (l + X± )( ] "H r ) 

J^cyl c:yl J "cyt 

2 Xa [ n r n z m H r ' m H z dv + [ \( ex H z ) 2 + ( ex H r ) 2 } dv - [ H 2 dv 
Jn cyl JR3\n cy i 1 J Jr3 

[ (H + m H l z fdv - [ H 2 dv + f \(H + m Hl) 2 + (1 + x ±) ( "H})' 

a [ n r n z (H + '"Hi) ' m H}dv + [ \(H + ex Hl) 2 + ( ex H^,) 2 } dv - [ H 2 c 

Jn cyl J^ 3 \n cyl 1 J Jr 3 

j {(1 + X,|) [2tf ' m H] + (' n Hl) 2 ] + (1 + x± ) ( m H}.) 2 } dv + 



dv + 
2 dv 



- X\\Hq [ I dr 

cyl 



dv + 



cyl 



2 Xa H ! n r n z ' m H l r dv + [ {2H ex i^ + ^H l z f + f x H l r f\ 
Jn cyl JR 3 \n cyl J 



dv 



(1 + X||) / \2H m Hl + ^Hl) 2 ]dv + {l + Xi _) / (>"H 1 r ) 2 dv + 



cyl 



cyl 



2 Xa H I n r n z in tf> + / \2H ex i^ + ( ex ^) 2 + {"H}) 2 } 
Jn cyl JR 3 \n cyl 1 J 



dv 



= I {(l + Xll )C"H 1 z ) 2 + (l + x± )(>"Hl) 2 }dv+ f {(**Hl) 2 + (**Hl) 2 }dv + 



2H I Xa I n r n z ' m H}dv + (1 + X ||) / m H 



cyl 



cyl 



Jm 



ex Hldv 



For the aims, discussed in section 3.1, we also give the linear in Xa representations 
U = U + Xa U 1 , ]n Hl=^H}.+^l, ' m H l z =^Hl + ^l 

where 



(A2) 



Uo = (1 + X±) 



Ui 



2fl- ((l + X±) / ^Hldv+ [ "Hldv) , 

V J^cyl J^ 3 \Vcyl J 

2 ^ + X± ^ [ (^H~l '™H\ + '"fl? + #o ] "Hl) dv + 



( ex ^) Z + ( ex i? r 1 ) 2 ' du+ (A3) 



(A4) 



cyl 
2 



+ 2n r n 2 tf in #i + 2H ™H\ 



In the last formulas we introduced the following notations 



Ad 



^H} = -kAxhtkr) sin kz , m H} = - Xa -^-N 2 h{kr)smkz , 

A d 



(A5) 



"H\ = -kA^oikr) cos kz , m H l z = - Xa -rNiI (kr) cos kz 

k 
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In fact, Uq was calculated in [4] 

U = -ttL ( X ±(okRH ) 2 I (kR)K (kR) . (A6) 
Calculate the integrals in (A4) taking in mind x± *C 1 

— f ^H\^H\dv = AxAfN-L ! Il{kr) cos 2 kzdv='^A 1 A 4> N 1 R 2 [l^(kR)-lf(kR)] 

= -^xdCoH^kRfN, [l 2 {kR) - I 2 {kR)] , (A7) 

— I ^H}'""Hldv = A 1 A 4> N 2 [ l\ (kr) sin 2 kzdv = ^A l A (j> N 2 R 2 [I 2 (kR) - I (kR)I 2 (kR)] 

Xa JQ cyl J^lcyl 

= -^XL(C H ) 2 (kRfN 2 [ll{kR) - h{kR)I 2 {kR)] , (A8) 

— / ' m H]dv = — r^Ho I (kr)coskzdv = TrLCoR-^N 1 H Io(kR) = 

Xa Jn cyl k Jn cyl k 



^2 ; n Io(kR) 



-nLmiCoHoYkRj^ , (A9) 



f ( in i?i) 2 dv = k 2 A\ j Ilikr) cos 2 kzdv = ^A 2 {kR) 2 [^{kR) - I 2 (kR)] 
ttL 

= —(x±CoH ) 2 (kR) A [l 2 (kR)-I 2 (kR)]K 2 (kR), (A10) 

H / n r n z m Hldv = -k 2 A X A 2 H Q I I 2 (kr) sin 2 kzdv 

J Q, C yi J Qcyl 

= -^A 1 A 2 (kR) 2 H [if(kR) - I (kR)I 2 (kR)] 

= ^X^Ho) 2 {kRf[l 2 {kR)-h{kR)I 2 {kR)}^^, (All) 

H ' m H}dv = -kAiH G I I {kr) cos kzdv = TvLAiH C kRI (kR) 

= irLx±((oHo) 2 (kR) 2 I (kR)K (kR) . (A12) 



It is quite surprising that among all above integrals there is only one (A9) which dominates over the 
others in the region x± ^ 1- Thus, in this limit we finally have 



2 , D Io(kR) ( Dr . . . . nAT Xa 



-U = 7r L(CoH rx±kR^—^ X ±kRh(kR)K (kR) + 2^— • (A13) 

hikR) V X±J 

The last expression shows that the anisotropic part of U can prevail over the isotropic one provided 
that, 

\NiXa\ > xl • (A14) 
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B On the solution of Equation (43). 



The non-homogeneous equation 

^ + * = CI Q {x) , *(x) = CG(x)I (x), (x ) = 

leads to the non-homogeneous equation for the amplitude function G(x) 

d 2 G ( h(x) 1\ dG , 

+ 2 -F7i( + - hir = 1 . G(x o ) = 0, 



cfa; 2 \ /o(^) x J dx 
or, after substitution dG/dx = S(x), essentially simplifies the problem 

ax \ 1q{x) x J 

Its solution reads 



«*> - -(-/( 2 lKHM/(^ 

and finally 

i ril{ y ) 



m + l -\dt 



(t) 



(y) 



ydy 
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